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In the current work spectroscopy and possibility of observation at the LHC of tetraquarks com- 
posed of 4 heavy quarks is discussed. Tetraquarks concerned are — [cc][cc], T±b = [66] [55] and 
T2[bc] = [&c][bc]. By solving nonrelativistic Schroedinger equation masses of these states are found 
with the hyperfine splitting accounted for. It is shown that masses of tensor tetraquarks T4 C (2 ++ ) 
and T 2 [i)c](2 ++ ) are high enough to observe these states as peaks in the invariant mass distribu- 
tions of heavy quarkonia pairs in pp — > + X — > 2J/ip + X, pp — > T 2 [b c ] + X — >• 2B C + X and 
pp — i> T 2 [b c ] + X —¥ J/4>T(1S) + X channels while Tib is under the threshold of decay into a vector 
bottomonia pair. 
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O , I- INTRODUCTION 

, Recent observation of J/?/>-meson pairs production in proton-proton collisions at 7 TeV energy at LHC renewed 
' interest to the 4 heavy quarks final states. In the low invariant mass region these quarks can form bound states 
(called tetraquarks) which can be produced in hadronic experiments. Thereby we would like to discuss physics of 
these states, elaborate their mass spectrum and possibility of experimental observation. 

Conception of tetraquarks, i.e. mesons composed of 4 valent quarks (qqqq), was first introduced in works 0, [IJ m 
II ■ 1976. For example, ao-meson and cr-mesons was treated as possible tetraquark candidate However, it is hard 

Oh' to determine quark composition of a particle in the light meson domain so these ideas were not developed further. 
Observation of new unexpected states such as A (3872) 0, § gave this idea a new impetus [H, Eccentricity of 
these particles consists in fact that according to the modes of their production and decay they contain a cc pair but 
they can not be included in the well known systematics of charmonia. Later similar particles were also found in the 
C^l 1 bottomonia sector [ldl - fl4l |. It is natural to ascribe these mesons to tetraquarks (QqQq), where q are light and Q - 
heavy (6 or c) quarks. 

However, situation when all quarks composing a tetraquark are heavy was not treated in details yet. This possi- 
bility seems to be quite interesting as in this case determination of meson's quark composition becomes simpler and 
its parameters can be determined by solving nonrelativistic Schroedinger equation. Our work is devoted to these 
particular questions. 

In our recent paper [15] we consider tetraquark T± c = [cc\[cc\ in the framework of diquark model. The hyperfine 
splitting in that paper was described through interaction of total diquark spins. Now we would like to study also 
tetraquarks = [bb][bb] and T 2 [bc] = [6c] [6c]. The last state is especially interesting since, in contrast to tetraquarks 
build from four identical quarks, both singlet and triplet spin states of the diquark are possible. It is clear, that 
hyperfine interaction of spin-singlet diquark cannot be described with the method used in our previous paper, so some 
other approach should be applied. 

In the following section spectroscopy of (cccc), (6666) and (bcbc) tetraquarks is concerned. Possibility of observation 
of these particles in hadronic experiments is discussed in the third section. We summarize our results in the short 
conclusion. 

II. SPECTROSCOPY 

A. General preliminaries 

In the current study diquark model of tetraquark is used. According to this approach tetraquark 

T = Q 1 Q 2 Q 3 Q i 
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consists of 2 almost point-like diquarks T>i2 — [Q1Q2] and P34 = [Q3Q4] with certain quantum numbers (such as 
angular momentum, spin, color) and mass. What concerns color configuration, two quarks in diquark can be in 
triplet or sextet color state. According to the manuscript [IB] in the sextet configuration diquarks experience mutual 
repulsion so we restrict ourselves to the (anti)triplet color configurations. Angular momentum of the diquark system 
equals in the ground state, so its spin is equal to the sum of quark spins which is or 1. What concerns diquark 
mass, it can be determined by solving Schroedinger equation with a correctly selected potential. According to work 
[l7T ] it is possible to use quark-antiquark interaction potential used in heavy quarkonia calculations with additional 
factor 1/2 due to the different color structures. 

In the diquark model tetraquark mass can be determined by solving 2-particle Schroedinger equation with point- 
like diquarks. As [Q1Q2] and [Q3Q4] diquarks are in (anti)triplet color configuration, potential of their interaction 
coincides with that of quark and antiquark in heavy quarkonia. Hyperfine splitting in this system can be described 
by the hamiltonian Q 

= M + 2^, % (S i S J -), (1) 

i<j 

where Mo is the tetraquark mass without splitting, Si — spin operator of i-th (anti)quark, and Kij are constants 
determined from experimental data analysis or theoretically. When dealing with potential models, Kij coefficient can 
be obtained from the value of the QiQj system wave function at origin: 

when both Qi and Qj are quarks or antiquarks in the color triplet state, and 

' * a s \R (ij) (0)| 2 (3) 



2 9m Ql m Q2 



if Qi and Qj are quark and antiquark in color singlet configuration. Three loop expression for the strong coupling 
constant was used when calculating Kij constants [17] and scale for it was taken equal 

^2 = 2m Ql m Q2 
m<3i + m Q2 

where (T) is average kinetic energy of quarks which is equal 

(T d ) = 0.19 GeV 

and 

(T„) = 0.38 GeV 

for the triplet and singlet states respectively. Values of diquark wave functions in origin are presented in manuscript 
[l7| while for mesons they can be calculated from the leptonic width T ee or leptonic constant / of meson in question: 



\R(0)\ 



2 _ J_ r _ M 2 f 
a 2 ef ee 4 ~ 9 ' 



Up to the hyperfine splitting tetraquark can be described by its total spin J, spins of diquarks S12 and 5*34 
constituting it and its spatial and charge parities P and C: 



|0++) = |0;0,0), |0++') = |0;1,1) 



1+- ) = 11:1,1) 



|1 +± ) = ^(|1;0,1)±|1;0,1», 
|2++) = |2;1,1> 

In this treatment all states are confluent with mass Mq. If spin-spin interaction is accounted for, masses of |1 ++ ) and 
1 2++) states shift: 

M (1++) = (1++ \R\ 1++) = M - k 12 - k_, 
M (2++) = 2m [12] + K12 + k + , 
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where following designations are introduced: 



K± 



2k 14 ± « 13 ± K 24 



and |0 ++ ), |0 ++ ) and |1 + ), |1 + ) states mix with each other. In the scalar tetraquarks case mixing matrix has the 
following form: 



H 



r |o++> - 






0++') 





Mo — 3«i2 — V3k- 

— \/3k_ Mq + Kl2 — 2«_| 



|0++> 
0+ 



and for |1 + ) tetraquarks: 



H 



\ ll + -> " 









Mq - K\2 + K13 - K 2 4 

«13 - ^24 Mq + K 12 ~ K+ 



+ - 
11+-' 



B. [QQ][QQ] 

In the case when quarks of the same flavor are involved, Fermi-Dirac statistics leads to the additional restrictions 
on the diquark quantum numbers. Indeed, permutation of quark indices should change sign of the total diquark wave 
function. As quarks are in the antitriplet color state, color part of this function is antisymmetric. Radial wave function 
is symmetric as quarks are in the S'-wave. So spin part of the wave function is to be symmetric too. Consequently 
the total spin of the S'-wave diquark can only be equal to 1. As a result only |0 ++ ), |1 ++ ) and |2 ++ ) diquark states 
remain. They do not mix with each other after the spin-spin interaction is accounted for. Masses of these states are 
equal 

M (0 ++ ') = M + k 12 - 2k+, 
M(l+-') = Mq + Ki2 - K+, 

M (2+ + ) = Mq + K 12 + K + . 

It worth mentioning that this splitting scheme agrees with the result of the work [15] in which interaction of the total 
diquark spins was concerned. 

To obtain numerical values of tetraquark masses one needs to know the unsplitted mass Mq and coefficients Kij in 
the Hamiltonian (JTJ. These coefficients can be calculated using expressions © and (|3J). The following values of quark 
masses were used: 

m c = 1.468 GeV, m b = 4.873 GeV. 

Diquark masses without hyperfine splitting are given in (l7l | while Mq mass of the tetraquark was calculated using a 
procedure similar to that described in [l7j . 

Let us begin with the tetraquark composed of 4 c-quarks, T 4c — [cc\[cc\. Mass of a ground state and value of radial 
wave function at origin for a [cc] diquark given in [17| are 

m [cc] = 3.13 GeV, i? [cc] (0) = 0.523 GeV 3/2 . 
Value of radial wave function at origin of the (cc) state determined from the leptonic width of J / ip- meson equals 

#(cc)(0) = 0.75GeV 3/2 . 
Spin-spin interaction coefficients calculated using expressions ([2]) and ([3]) are equal 

«12 = ^34 = K[ C c] = 12.8 MeV, 

«13 = K 23 = «i4 = K 24 = «(cc) = 42.8 MeV. (4) 
Without hyperfine splitting T 4c tetraquark mass equals 

M = 6.124 GeV, 
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and with it this state splits into scalar, axial and tensor mesons with masses 



0++ : 


M = 


5.966 GeV, 


M 


-Afth 


= -228. MeV, 




M = 


6.051 GeV, 


M 


- Ma 


= -142. MeV, 


2++ : 


M = 


6.223 GeV, 


M 


-Ma, 


= 29.5 MeV. 



In expressions above differences between the tetraquark masses and a J/ip-meson pair formation threshold are also 
noted. It can be seen that only tensor state lies above this threshold and can be observed in the T4 C (2 ++ ) — > 2J/tp 
mode. It worth mentioning that scalar tetraquark which is slightly under the J/-0-pair threshold can however decay 
by the T4 C (0 ++ ) — > (J / %/))* J / ip — > n + n~J/tp channel i.e. with one J / ip-meson being virtual. So it can be observed 
as a peak in the J/ip invariant mass distribution. 

For a tetraquark built from 4 6-quarks, i.e. T± b — [66] [66], situation is entirely similar to the previous case. Mass of 
the [bb] diquark and values of radial wave function at origin for it and for the (66) ground state are 

m [bb] = 9.72 GeV, % fc ](0) = 1.35GeV 3/2 , %5)(0) = 2.27GeV 3/2 . 
Mass of the T^ b tetraquark without hyperfine splitting is equal to 

M = 18.857 GeV, 

and spin-spin interaction coefficients are 

«12 = k 34 = K [bb] = 5.52 MeV, 

«13 = K 2 3 = KU = «(65) = 27 - 1 MeV ' ( 5 ) 

With hyperfine splitting one obtains the following masses of the T^ b states: 

0++' : M = 18.754 GeV, M - M th = -544. MeV, 
1+"' : M = 18.808 GeV, M - M th = -490. MeV, 
2++ : M = 18.916 GeV, M - M th = -382. MeV. 

It can be seen that in this case all the states are under the T(IS') pair production threshold M t h = 2m X (is)- 

C. [bc][bc] 

The situation is more interesting in the T 2 [5 C ] = [6c] [6c] tetraquark case. In this case [6c] diquark spin can be or 1 
and all states mentioned in Section II-A exist. Diquark mass and value of its radial wave function at origin are |17| 

m [bc] = 6.45 GeV R [bc] (0) = 0.722 GeV 3/2 . 
Radial wave function at origin for the color singlet (6c) state can be determined by its leptonic constant fg — 500 MeV 

m- 

R m = 1.29GeV 3/2 . 

So spin-spin interaction coefficients are equal 

K12 = k 34 = K[ bc j = 6.43 MeV, 
«14 = k 2 3 = K(be) = 27.1 MeV. 

Values of K13 = K/ bb \ and K24 = K( cE ) constants were given in expressions Q and ([5]). Without hyperfine splitting 
T 2 [ bc ] tetraquark mass equals 

M = 12.491 GeV, 

and with spin-spin interaction accounted for this state splits (see Fig[T]) into 
• Two scalar states with masses 

++ a : M = 12.359 GeV, M - M th = -191. MeV 
0++6 : M = 12.471 GeV, M - M th = -78.7 MeV, 
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Figure 1. [&c][6c] tetraquark mass spectrum 
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• Two 1 + states with masses 

l+"a : M = 12.424 GeV, M - M th = -126. MeV 

: M = 12.488 GeV, M - M th = -62.5 MeV, 

• One 1 ++ meson with mass 

1++ : M = 12.485 GeV, M - M th = -64.9 MeV, 

• One tensor meson with mass 

2++ : M = 12.566 GeV, M - M th = 16.1 MeV. 



Mass of the two B c mesons is selected for the threshold value in these expressions, M t h = 2ms c = 12.55 GeV. It can 
be seen that only tensor tetraquark T 2 [bc](2 ++ ) lies above this threshold and thus can be observed as a peak in the 
_B c -meson pair invariant mass distribution. 

In paper [l8[ tetraquark states were also considered in the framework of diquark model. Picture of hyperfine 
splittings of T 2 [b c ] -tetraquark, presented in this paper is in good agreement with our results. Predictions for masses, 
on the other hand, are about 700 MeV higher, than our values. As a result, according to this paper all tetraqurk 
states should lie above 2B* and J/ipT thresholds. We think, that the main reason for difference between these two 
works is the neglection of binding energy in tetraquark tetraquark and diquark spectra, that is negative. For example, 
for tetraquark state before hyperfine splitting we have SE — Mq — 2m\p c \ k, 410 MeV. 



III. PRODUCTION 
A. Duality relations 



Duality relations can be used to estimate production cross sections of the particles in question. Let us consider 
formation of two diquarks in gluon interaction gg — > [Qi^HQsQ^- Above the 2 doubly-heavy baryon production 
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Figure 2. Invariant mass distribution of J/ip-meson pairs with expected Ta c tetraquark contribution. Dashed vertical line 
corresponds to the two doubly-heavy baryons formation threshold 2Mh cc . 



threshold these diquarks can hadronize into 4 open heavy flavor mesons, 2 heavy quarkonia or form a bound state, 
i.e. tetraquark. According to our estimations this tetraquark would preferably decay into a vector quarkonia pair. 
Indeed, decay into light mesons is suppressed by the Zweig rule, production of 4 open heavy flavor mesons is prohibited 
kinematically and formation of pseudoscalar quarkonia requires flip of the heavy quark spin. That is why the following 
duality relation can be written: 

2M *QQ 2M *QQ 

S T = J dm gg a\gg->T->2Q]=e J dm gg a(gg -> [QiQ 2 ] + [Q 3 Qa]), (6) 

2Mq 2m [QQ] 

where e factor stands for the other possible decay modes. This value is to be compared with the integrated non- 
resonant cross section of quarkonia pairs production in the same duality window: 



S 2 q = J dm gg o [gg -> 2Q] . (7) 

2A/ e 

As tetraquark states are typically narrow, these mesons can be observed as peaks in the quarkonia pairs invariant 
mass distributions despite that St *C S%q relation holds. As tetraquark width is small compared to the detector 
resolution A w 50MeV, its peak can be modeled by a Gaussian form with corresponding width. Therefore cross 
section of the gg — > T — > 2Q Breit-Wigner process is replaced with the following expression: 

where preexponential factor is selected according to the duality relation {BJ. 



B. T 4Q 

Let us begin with the tetraquarks composed of the 4 identical quarks. In the T^ c case only tensor state lies above 
the 2 vector charmonia production threshold. Integrated cross sections calculated using expressions © and J7]) arc 
equal 

5 Tlc = 0.7 pb GeV S 2J/i , = 20 pb GeV, 

where suppression factor is selected to be e = 0.2. Invariant mass distribution for the J/^-meson pairs with expected 
T^c tetraquark contribution is shown in Fig. [2] 

As already mentioned, in the tetraquark case even tensor state is under the two vector bottomonia formation 
threshold so its observation in their invariant mass distribution is doubtful. 
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Figure 3. Invariant mass distribution of _B c -meson pairs (left plot) and J/ipT(lS) (right plot) with expected contribution of 
the T 2 [b c ] tetraquark. Dashed vertical line corresponds to the two doubly-heavy baryons formation threshold 2M~ i)e . 



C. T 2 [b c ] 

Let us turn to the T 2 \bc] = \bc\\bc\ tetraquark. In this case pseudoscalar B c -meson and vector i3*-meson decaying 
into B c "f can be experimentally observed. Thus tensor tetraquark T 2 [bc](2 ++ ) can be observed as a peak in the B c - 
meson pairs invariant mass distribution. However T 2 [f, c ](2 ++ ) — > 2B C decay requires flip of the heavy quark spin and 
is suppressed by the factor 

t i [££i ~ 2.6 x 1CT 3 . 

m[bc] 

Integrated cross sections ((6j) and (0 are equal 

S T2lba] = 0.13 fb GeV S 2Ba = 6 fb GeV, 

where M~, bc =6.82 GeV [l{| is used. Invariant mass distribution of the £? c -meson pairs with expected contribution of 
the T 2 [h c \ tetraquark is shown in Fig. [3^,. All T 2 [t c ] mesons lie under the B* pair production threshold. 

Decay of the T 2 ^ c ] tetraquark into the J/ipT(lS) vector quarkonia pair is also possible. In the color singlet model 
reaction gg — > J/ipT(lS) is prohibited so accounting for octet components of vector quarkonia is needed. This 
process was elaborated in the work [2(| results of which are used as the background for the tetraquark contribution. 
Suppression factor e~3x 10 -2 was used to obtain integrated cross sections (j6]) and (J7J): 

S T2[bc] = 1.5 fb GeV, S^r = 33fbGeV. 

Invariant mass distribution of the J / tpT (IS) pairs with expected contribution of the T 2 [h c ] tetraquark is shown in Fig. 
[3b. 



IV. CONCLUSION 



In our work tetraquarks composed of 4 heavy quarks are concerned. Spectroscopy of T^ c — [cc] [cc] , T^h = [bb] [66] 
and T 2 b c = [6c] [6c] states is elaborated and possibility of observation at the LHC is studied. 

Model used implies diquark structure of tetraquarks i.e. Q1Q2Q3Q4 tetraquark is assumed to consist of 2 almost 
point-like diquarks [Q1Q2] and [Q3Q4] being in triplet color configuration. In such a treatment tetraquark is entirely 
analogous to a doubly-heavy meson. So non-relativistic Schroedinger equation which gives reliable results for the 
charmonia and bottomonia description can be used to obtain its parameters. 

In the case of T^ c and T^h tetraquarks Fermi principle imposes constraints on the possible quantum numbers of 
diquarks: in the antitriplet color configuration total spin of the .5- wave diquark can be only equal to 1. So with the 
hyperfine splitting accounted for 3 states with the following quantum numbers appear: ++ , l" 1 and 2 ++ . In the 
charmed tetraquark case only tensor meson lies above the vector meson pair formation threshold and its peak can be 
observed at LHC in their invariant mass distribution. For the T^h tetraquark all the states lie under the T(IS') pair 
formation threshold. Both and 1 spins are possible for the diquark in the T 2 h c tetraquark. So 6 states arise after 
the hyperfine splitting of the unaffected state: two of ++ , one 1 ++ , two of l" 1 and one 2 ++ . 
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In the last section we estimated the possibility to observe tetraquarks concerned in the inclusive reactions gg — > 
T — > 2Q. According to our calculations T^ c {2 ++ ) tensor state can be observed as a peak in the J/^-meson pairs 
invariant mass distribution. T 2 [6 C ](2 ++ ) tetraquark can be observed in both 2B C and J/?pT(lS) modes. 
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